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Low energy excitation of surface states of a three-dimensional topological insulator (3DTI) can be 
described by Dirac fermions. By using a tight-binding model, the transport properties of the surface 
states in a uniform magnetic field is investigated. It is found that chiral surface states parallel to the 

2 

magnetic field are responsible to the quantized Hall (QH) conductance (2n+ 1) V multiplied by the 
number of Dirac cones. Due to the two-dimension (2D) nature of the surface states, the robustness 
of the QH conductance against impurity scattering is determined by the oddness and evenness of 
the Dirac cone number. An experimental setup for transport measurement is proposed. 

PACS numbers: 71.70.Di, 72.10.-d, 73.20. At, 73.43.-f 



I. INTRODUCTION 

Topological insulator (TI), a new quantum state in 
which the bulk is an insulator and the surface is metallic, 
has attracted intensive attention^ since its theoretical 
predictions 0-|3| and its experimental realizations in 
three dimensions. Unlike normal insulators, the bound- 
ary (edges or surfaces) of a finite TI sample supports 
extended and current carrying states in its bulk gap. In 
the presence of time reversal symmetry, TIs can be clas- 
sified by a Z2 invariant integer i/[9l4llj into a weak TI 
for even number of the Dirac cones (called v = 1) and 
a strong TI for odd number of the Dirac cones (called 
v = —1). The edge/surface states of a weak (strong) TI 
can (cannot) be destroyed by defect scatterings d, y|. 

The surface states of a three-dimensional topologi- 
cal insulator (3DTI) are described by two-dimensional 
(2D) Dirac cones centered at the time-reversal invariance 
points in the Brillouin zone. 2D Dirac fermions with a 
single cone are predicted to have many novel properties 
such as the Klein paradox, anti- localization, etc. [12i4l5l | 
In the quantum Hall effect (QHE) regime (strong mag- 
netic field), the Landau levels of a single cone Dirac 
fcrmion is unevenly distributed 



n = 0,±1,±2,. 



(1) 



This distribution has been confirmed by several 
experiments 16-ljj. Another interesting prediction of 
one Dirac cone fermion is the half (in the unit of e 2 /h) 
quantized Hall (QH) conductance 14, lj| 2(| because one 
Dirac cone carries a it Berry phase. Although many prop- 
erties of strong 3DTIs have been confirmed since its ex- 
perimental realization in several semiconductors, the half 
QH conductances has not been observed in experiments 



yet, and this will be our focus. 

In this paper, we investigate a 3DTI in a uniform mag- 
netic field within a well-known tight binding model. Our 
calculations yield indeed the well-known Landau levels 
(Eq. (JT])). However, the issue of the half QH conductance 
is subtle. All surface states of a finite 3DTI are fully con- 
nected with each other @, d, 21 1 and surface states living 
on the sample surface may move from one side to an- 
other. In other words, the well-defined Hall voltage mea- 
surement in usual 2D electron gases (2DEG) is not so 
clear in the 3DTI cases since two surfaces cannot be sep- 
arated by the bulk of a finite 3DTI as it did in 2D QHE 
systems. The current does not come from ID channels, 
but from 2D surface states. The structures of these side 
surface states are more complicated than the edge states 
in 2DEG, therefore well-defined quantum Hall plateaus 
(2n + 1)tt can only be defined from transverse current 
instead of the Hall voltage. Our tight-binding calcula- 
tions confirm the earlier predictions based on an effective 
theory of Dirac fermion in curved 2D spaces [2lJ- The ro- 
bustness of the Hall conductance against impurities is 
also studied. Our tight-binding model and their basic 
properties are introduced in the next section. Numerical 
results and their discussions are presented in section III, 
followed by the conclusion section. 



II. MODEL 



A well-known tight-binding model of 3DTI[22[ is 



H (k) = e (k)I ix4 + J2d a (k)T a 

a=l 

d a (k) = (A 2 sink x ,A 2 sink y , A 1 sink z ,M(k),0), (2) 



2 



where e (k) = C + 2Lq + AD 2 - 2D X cos k z - 2D 2 (cos k x + 
cos k y ),M = M - 2Bi - 4P 2 + 2B 1 cos k z + 2P 2 (cos k x + 
cosky) and r 1>2 ' 3 ' 4 ' 5 = (a x ® s x , a x ® s y ,a y ® I 2X 2, (J z ® 
^2x2, <Ti®s 2 ) in the basis of four states (|P1+, f), |P1^ ,i 
), |P2~,f), |P2~,J.)). After inverse Fourier transforma- 
tion of equation the real space version of this model 
on a cubic lattice can be written in the form 



(ij)ss 



tij c ls C j 



H.c. 



(3) 



where i,j are lattice site indices, s € {|P1+, t), |Plj~,4- 
), |P2j,t), |P2~,4)} is the spin-orbital index. The first 
term in equation is the on-site energy, the second 
and third terms describe the hoppings between nearest 
neighbor sites. The effect of non-magnetic impurities can 
be included by adding a term 



Hr 



(4) 



where Vi S distributes randomly in the energy range 
of (-W/2,W/2). The magnetic field B is intro- 
duced thr oug h the Peierls substitution of hopping 
coefficients [23, [2J] 



exp( 



2iri 



dl ■ A)t 



(5) 



where 0o = e/h. 

In the clean limit (H] = 0), this 3D model is fully 
gapped in the energy range (— M, M). The proper- 
ties of surface states within this gap are determined by 
the Z 2 topological numbers related to the time reversal 
polarizations [3, 0, 0] 



where vf = A 2 \j\ — (-§x) 2 and h = 1 is adopted. 

Before presenting calculations of tight binding models, 
let us first look at effective model (Eq. 0) in a field 
B = (0, 0, B) perpendicular to the surface. In the Landau 
gauge A = (-By, 0,0), H xy is[H 



H X y Vp[<J X ( K ky Ay) <Jy{k X -^c)] 

k y + i(k x + eBy) 



= v F 



i(k x + eBy) 







, (9) 



where fcj = — idi. The eigenvalues are Landau levels [13| 

E xy = ±v F V2neB, n = 0,1,2,- •• (10) 

However, for surface parallel to magnetic field, say, x — 
z plane, the Hamiltonian is 

H xz = a x (p z - A z ) - a y (p x - A x ) (11) 
p z + i(p x + eBy) 



Pz -i(Px + eBy) 







(12 



Since eBz commutes with k x and k yi the only difference 
between Eq. (fT2| and Eq. is a shift —eBy of the Dirac 
point in k x direction. The eigenvalue shows a shifted 
Dirac cone 



E xz = ±v F \/{k x + eBy) 2 + k 2 z . 



(13) 



We will see that states in the surfaces parallel to the mag- 
netic field play an important role in the QHE in 3DTI. 



III. NUMERICAL RESULTS 



y/detHr,)] 
pf[ w (ro] 



-sgn(d 4 (k = r<)) = ±1, (6) 



at 8 time-reversal invariant points in the first Brillouin 
zone: ^,3,3,4,5,6,7,8 = (0,0,0), fa 0,0), (0,tt,0), (0,0, tt), 
(0,7r,7r), fa0,7r), fa7r,0), fa7r,7r), where (k) = 

(u-k,m|0|uk,n)< K v = ]li=i = tri e system is a 
strong 3DTI with odd numbers of the Dirac cones on each 
surface. Otherwise, the system (v = 1) is a weak 3DTI 
with even numbers of the Dirac cones on each surface. 
In terms of model parameters, the system is in a strong 
(weak) TI phase when B % /M > 1 (P^/M < -1)@. 

The surface states around each Dirac cone is approx- 
imately described by the Dirac fermion Hamiltonian (z- 
axis is normal to the surface) 0, El HE Hi 



H xy = VF{cr x k y - 
with a linear dispersion relation 



<jyk x ) 



(7) 



E xy = ±v F \ k% + k 2 



To be specific, we consider Hamiltonian j3| on a cubic 
lattice with size N x x N y x N z . For a slab 00 x 00 x N z 
of thickness N z in z-direction and infinite in other two 
directions 0, the surface states of the TI can be dis- 
played by plotting the dispersion relation E(k x ,k v ) in- 
side the bulk gap so that this dispersion relation must 
be from the surface states. Fig. Q] (a) is the disper- 
sion relation of a slab without a magnetic field. The 
Dirac cone can be clearly seen. There are two degen- 
erated Dirac cones located on lower (z = 1) and up- 
per (z = N z ) surfaces, respectively. The existence of 
such surface states roots on the topological property of 
time reversal symmetric systems [H The magnetic 
field breaks time reversal symmetry and the fate of sur- 
face states will be tested in the following tight-binding 
calculation. A uniform magnetic field generally also 
breaks the lattice translational symmetry 26[ and results 
in the Hofstadter butterfly spectrum [27|. If the mag- 
netic flux through a unit cell is a fraction of the flux 
quanta, Ba 2 — ^(f>o, where p and q are two prime num- 
(8) bers and a is the lattice constant, periodic structure is 
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FIG. 1: (Color online) Dissipation relations of a strong 3DTI 
with model parameters A\ = A% = 1, B\ = B<z = 1, C = 0, 
Di = D 2 = and M = 0.4. (a): E(k x ,k y = 0) for a slab with 
geometry oo x oo x 40 and without magnetic field, (b): Similar 
to (a), but in a perpendicular magnetic field B z = jf^ a i ■ Inset 
of (b): the Landau levels E n versus \fn. (c): E(k x ,k z = 0) 
for a slab with geometry oo x 60 x oo in a parallel magnetic 
field B z = (d) E(k x ) for a bar of oo x 60 x 40 with 



periodic boundary condition in z direction, in a magnetic field 
B z — 4 QQ a i , in contrast to its infinite-size version (c). 



restored with an enlarged and y-elongated unit cell of q- 
times of the original one, a' x = a x and a' y = qa y , then 
k x € [— 7r/a,ir/a and k v G [— ir/(qa), — 7r/(ga)] are good 
quantum numbers [28l |29||. Fig. Q] (b) is the dispersion 
relation of the same slab as that for Fig. [1] (a) in a per- 
pendicular magnetic field B = (0, 0, B). Instead of linear 
dispersion relation in zero field, discrete Landau levels 
E n appear. We have confirmed that there are an upper- 
lower surface double degeneracy and all wavefunctions 
are strongly confined inside (top/bottom) surfaces. Al- 
though the time-reversal symmetry is broken, the mag- 
netic field does not destroy the 2D nature of the sur- 
face states. The dependence E n ~ y/n is also verified as 
shown in the inset of Fig. Q](b). Such Landau levels have 
already been observed in recent experiments 171 118|. 

We plot also the dispersion relation for the infinite x—z 
surfaces in a parallel magnetic field B = (0, 0, B) in Fig. 
Q] (c) . Shifted Dirac cones corresponding to two surfaces 
(y = 1 and y = N y ) can be clearly seen, as predicted in 
Eq. (|T3|) . We will see that these side surface states play 
an important role in the following discussion. Fig. [T] (d) 
is similar to (c) when the infinite surface is confined along 
z-direction with a periodic boundary condition that elim- 
inates possible edge states. Subband structures appear 
seen due to the z-confinement. 

With quantized Landau levels, one naturally expects 
quantum Hall effects. In conventional 2DEGs, quan- 
tized Hall conductance comes from the spatially sepa- 
rated counter-propagating edge channel(s) on the two 



sides of a samples [301]. The situation in 3D is subtle be- 
cause a surface state may cover all surfaces because they 
are fully connected @, d, [2l|. In other words, a surface 
state cannot be confined on only one surface, resulting in 
more interesting with richer physics. 

To study possible current carrying surface states, we 
consider a bar geometry of finite widths of N y and N z 
in both y— and z— directions. The dispersion relation 
along x— direction in a magnetic field is plotted in Fig. [5] 
The constant energy in the middle of the curves shows 
that the Landau levels exist in the middle of the sur- 
face. The Landau levels float up near the edges (junc- 
tions of two surfaces). We have confirmed that the left- 
most (k x < 0.1) and right-most (k x > 0.8) parts of the 
curves correspond to states distributed mostly on the left 
(y = 1) and the right (y = N y ) side surfaces, respectively 
(Fig. [3]) . These side surface states are spin-polarized due 
to spin-orbit coupling |21|. as plotted in Fig. |4] 

As illustrated in Fig. [HJa) to (e), varying k x from ei- 
ther side of the curves to the middle flat parts, the states 
move from the side surfaces of y = 1 ant y = N y (Fig. Q] 
(d)), to the top and bottom surfaces of z = 1 and z = N z 
(Fig. H(b)) continuously. This is possible because a par- 
ticle can move from the side surfaces to the top/bottom 
ones without passing through the bulk. In other words, 
the surface states live essentially on a closed surface of a 
finite 3DTI sample. One cannot have ID edge channels 
as in the 2D case. The discreteness of these states simply 
originates from the finite thickness N z in z-direction, as 
indicated in Fig. [2] (a) and (b). Increasing N z does not 
affect the Landau levels, but generates more side surface- 
state subbands. In the limit of large N z , these subbands 
pack densely and form two continuum cones, similar to 
the shifted Dirac cones in parallel magnetic field in Fig. 
□ (c). 

Similar to edge states in conventional 2D systems, due 
to non-zero dispersion, the side surface states can trans- 
port electrons in the presence of voltage along x direc- 
tion. The effects of side surface states have been no- 
ticed in a recent transport measurement (3l| . However, 
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FIG. 3: (Color online) Spatial distributions of wavefunctions on the y — z section of the bar for the states marked as red dots 
in Fig. [21(a). The magnetic field is in z-direction. (a) to (e) corresponds to state points a to e, respectively. 
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FIG. 5: (Color online) Schematic drawing (seen from above) 
of the active side surface channels at the Fermi energies Ei 
(a) and E2 (b) indicated by the orange lines in Fig. [2] (a). The 
channel numbers of forward moving rif and backward moving 
rib are indicated. 



FIG. 4: (Color online) Spatial distributions of local spin s z . 
(a) and (b) correspond to state points a and b in Fig. 
respectively. 



the side surface states are quite different from the edge 
states in 2DEG. For example, the dependence of E n (k x ) 
is non-monotonic on each side. This brings the coexis- 
tence of both forward and backward moving channels on 
each side, as illustrated in Fig. [SJ At each side and for a 
given energy, the numbers of forward moving channels Uf 



and backward moving depend on magnetic field B, as 
well as thickness N z of the sample. However, their differ- 
ence m — rib — 1) 3, 5, • • • is universal between any definite 
pair of adjacent Landau levels, as long as the dimension 
of the sample is not too small to couple states on any 
opposite surfaces [25|. A general theorem guarantees the 
existence of current carrying states 32[ . The directions 
of these net currents on two sides are reversed, thus they 
are chiral. Experimentally, such chiral currents can be 
measured by a multi-terminal measurement illustrated 
in Fig. |51 In such setups, the Landauer-Biittiker formula 
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FIG. 6: (Color online) The schematic drawing of 3D multi- 
terminal measurements of the TI (grey). The uniform mag- 
netic field B (red) is in the z— direction. The orange arrows 
represent the side surface channels that carry currents, (a) 
The 3D schematic diagram of a four-terminal measurement. 
A voltage V is applied to the left (1) and right (2) terminals. 
The potential of the front and back terminals are V/2 and 
V/2, respectively, (b) The vertical view of (a) from above, (c) 
The vertical view of a 3D six-terminal measurement, which is 
not preferred in our model. 



is Ij = J^i Tj^i [Vj — VA , where Tj^i is the transmission 
from lead i to lead j, Vi and It are the voltage and the 
current in the i— th lead respectively, satisfying Kirchoff 's 
law It = o[33j]. In the clean limit, all the active chan- 
nels are perfectly conducting and the transmission T}<_j 
can be read directly from the profile of the active sur- 
face channels, as shown in Fig. [B] For the four-terminal 
setup illustrated in Fig. [B] (a) and (b), straightforward 
calculations within this formalism give 
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(15) 



This G xy is the Hall conductance in the present problem. 
This is consistent with the previous calculation based on 
an effective model in curved space (2lj. but here under- 
stood within simple band theories. Similar results hold 



for the case of surface states with more than one Dirac 
cones and the only difference is the conductance is mul- 
tiplied by a factor of N conc , the number of Dirac cones. 
However, for a 3D six-terminal setup, shown in Fig. [5](c) 
and generalized from the traditional 2D Hall bar , similar 
calculations show that 
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which is generally sample dependent and not quantized. 
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FIG. 7: (Color online) Hall conductance G xy as a function of 
the Fermi energy E for various disorder strength W . G xy for a 
strong (a) and weak TI (b) is calculated from a four-terminal 
setup as illustrated in Fig. [5] The sample size is 38 x 18 and 
the model parameters are Ai = A2 = 1, B2 = 1, C = 0, 
Di = D 2 = 0, M = 0.8 and B x = ±1 for strong (weak) TI. 



One question arises: 
ductance G xy defined 



is the value of the Hall con- 
in Eq. (|T5|) robust against 
disorder 3J|? In zero magnetic field, it is well-known that 
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the surface states are robust for odd N conc (strong TI) 
and not for even N conc (weak TI) . Since charge transport 
is dominated by the side surface states at zero field, it is 
reasonable to expect that the robustness of G xy is also 
determined by the evenness and oddness of N conc . Specif- 
ically, for strong TI with single-cone surface states in a 
magnetic field, as illustrated in Fig. [2] forward and back- 
ward channels coexist at the same side, but they originate 
from a single Dirac cone and the backscattering is small. 
These are confirmed by numerical simulations of a four- 
terminal set-up as shown in Fig. [7] (a) for a strong 3DTI 
and in Fig. [7](b) for a weak 3DTI, by using the standard 
method of non-equilibrium Green's functions 33, 35| . In 
the weak TI case, the conductance G xy is twice as large 
due to double cone degeneracy. Also, its value is sensi- 
tive to disorder because of the inter-scattering between 
the two cones. Since a four-terminal numerical simula- 
tion is very resource consuming, in Fig. [7]we use a smaller 
size system. The physics remains the same as long as no 
direct coupling between surface states on two opposite 
surfaces. Quantized Hall plateaus predicted from above 
analysis, especially for the strong TI, can be clearly seen. 
In the presence of disorder, the quantum Hall plateaus of 
a strong TI (Fig. [7] (a)) survive much better than those 
of a weak TI (Fig. 0(b)). 



IV. SUMMARY 

Before ending this paper, it should be emphasized that 
we did not obtain the half conductance, or in a less strict 
sense, we saw the sum of half conductance from upper 
and lower surfaces (2 1|. since the Dirac fermion live on a 
closed and curved surface of 3DTI. By making an edge 
channel on one surface, one inevitably makes a new one 
on the opposite surface so that a Dirac fermion will not 
terminate somewhere on a surface. The exact half Hall 
conductance of one surface cannot directly be observed 
unless the surface state can be effectively confined in one 
isolated plane by some means. On the other hand, the 
quantized Hall conductance we obtained from the 2D sur- 
face states should not be viewed as a trivial application of 
Chern number theory of 2DEG[2^, [29[, since the surface 
states live on a closed 2D manifold embedded in a 3D 
space, which is topologically different from the 2D plane 
of 2DEG. 

In summary, we investigated the quantum Hall effect of 
a 3DTI in a magnetic field. The integer Hall conductance 
is carried by side surface states due to the non-separable 
nature of surface states enclosing the 3DTI. The quan- 
tum Hall conductance reflects the properties of side sur- 
face states that are parallel to the magnetic field. The 
quantum Hall effect thus offers a transport measurement 
to determine the topological property of a 3DTI: whether 
it is a weak or strong TI, the number of Dirac cones, etc. 
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